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ABSTRACT 

Aims. We investigate the relationships between three main optical galaxy observables (spectral properties, colours, and morphology), 
exploiting the data set provided by the COSMOS/zCOSMOS survey. The purpose of this paper is to define a simple galaxy classifi- 
cation cube, using a carefully selected sample of » 1000 galaxies. 

Methods. Using medium resolution spectra of the first Ik zCOSMOS-bright sample, optical photometry from the Subaru/COSMOS 
observations, and morphological measurements derived from ACS imaging, we analyze the properties of the galaxy population out to 
z ~ I. Applying three straightforward classification schemes (spectral, photometric, and morphological), we identify two main galaxy 
types, which appear to be linked to the bimodality of galaxy population. The three parametric classifications constitute the axes of a 
"classification cube". 

Results. A very good agreement exists between the classification from spectral data (quiescent/star-forming galaxies) and that based 
on colours (red/blue galaxies). The third parameter (morphology) is less well correlated with the first two: in fact a good correla- 
tion between the spectral classification and that based on morphological analysis (early-/late-type galaxies) is achieved only after 
partially complementing the morphological classification with additional colour information. Finally, analyzing the 3D-distribution 
of all galaxies in the sample, we find that about 85% of the galaxies show a fully concordant classification, being either quiescent, 
red, bulge-dominated galaxies (~ 20%) or star-forming, blue, disk-dominated galaxies (~ 65%). These results imply that the galaxy 
bimodality is a consistent behaviour both in morphology, colour and dominant stellar population, at least out to z ~ 1. 

Key words, galaxies: general - galaxies: evolution - galaxies: fundamental parameters (classification) 
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1. Introduction 

Galaxies have a wide variety of physical and observat ional p rop- 
erties. However, since the early works of Humason (1936) and 
Rubble ( 1938) it has been known that the morphologies of galax- 
ies correlate with their colours and consequently with the de- 
rived properties of their stellar population. Large efforts have 
been undertaken to establish reliable relationships between the 
morphologic al types, colours, and spec tral characteristic s of the 
galax ies (i.e. iFioc & Rocca-Volmerange, 1999; Bershadv et al.L 
I2OOOI) . The accurate determination of these relationships is of 
paramount importance to our understanding of galaxy forma- 
tion and evolution. The main galaxy optical properties are the 
spectral energy distributions, spectral line properties, internal 
velocity dispersions, morphologies, sizes, and structural com- 
ponents. From these observable quantities, we can derive in- 
trinsic physical properties, such as stellar and gaseous masses, 
metallicities, and star-formation histories; each of these as- 
pects provides crucial clues about how galaxies were created 
and have evolved. The most basic evolutionary processes and 
their corresponding observables are, however, not yet clear and 
we are still uncertain whether the correlations between galaxy 
properties observe d are driven by f undamental aspects of the 
galaxy population dConsehcel l2006h . The goal of this paper is 
to investigate the relationships between the three main optical 
galaxy observables (spectral properties, colours, and morphol- 
ogy) by taking advantage of the unique data set provided by the 
COSMOS/zCOSMOS survey. 

A fundamental property discovered by large surveys such 
as the Sloan Digital Sky Survey (SDSS) and the Two-Degree 
Field Galaxy Redshift Survey (2dFGRS), is that the optical 
colour distribution of galaxie s at low redshift can be repre- 
sented by a bimodal function dStrateva et al.L 120011: IWild et all 
l2005h . The bimodality separates the locus of star-forming galax- 
ies from that of early-type galaxies (ellipticals, SO, and early- 
type spirals), such that the fraction of the red, passive population 
increases progressively with the luminosity, becoming smaller 
than the blue fraction f or Mr > -20 and d ominating the dis- 
tribution for Mr < -21 dBaldrvet al.Ll2004 . The bimodahty is 
also present in different aspects of the galaxy population: struc- 
tural parameters such a s concentration index (Ball et al.. 2006, ) 
and disk-to-bulge ratio dDrorv & Fished | 2007h. and spect r al fea- 
tures such a s D4000 and Une ratios ('Kauff mann et aL ,_ 20031; 
iBalogh et aUl2004aHYan eTail l2006; Franzetti,etaLl |2007|), all 
exhibit strong bimodal distributions. The bimodality in galaxy 
propertie s, clearly evident in the local Universe, is still present 
at z ~ 1 dBell et al.L |2004|) . but it is less clear whether the two 
distinct galaxy populations hold at higher redsh i fts (up to z = ^ 
2; [G iallongo et al., 2005; C irasuolo et all 120071; ICassata et all 
I2OO8 :). A key test for galaxy evolution models is to explain why 
this bimodality occurs; it is still unclear which processes produce 
two sets of galaxies of different average colours and dispersions 
and different luminosity ranges. 

Determining the nature of the bimodality can be approached 
by means of large spectroscopic redshift surveys, studying the 
measured galaxy parameters as a function of redshift, luminos- 
ity, and spectral properties. The zCOSMOS survey opens the 
possibility to extend the study of bimodality out to z ~ 1 . 

zCOSMOS is a large redshift survey that is being undertaken 
in the COSMOS field as an ESO Large Programme (~ 600 hours 
of observation) with the VIMOS spectrograph at the VLT. The 
survey consists of two parts: the first part targets a magnitude- 
limited sample selected from ACS images to have Iab < 22.5 
and observed at medium resolution; this will provide spectra and 



redshifts for ^ 20, 000 galaxies out to z ~ 1.5 over two square 
degrees. The second and deeper part of the survey is acquiring 
spectroscopy to measure the redshifts of » 12,000 galaxies at 
1.4 < z < 2.5 at lower resolution and over the central square de- 
gree, by selecting targets using various colour-selection criteria. 

This paper, based mainly on the spectroscopic data of the 
zCOSMOS survey, takes advantage of the large potential of the 
COSMOS database. We present the analysis of a representative 
sample of the bright component of the full survey, consisting of 
around 1200 galaxies observed during the first year of obser- 
vations. Using the medium resolution spectra of these objects 
and complementing them with the optical p hotometry from the 
COSMOS multi-band catalog dCapak et al.i. .2007) and morpho- 
logical information derived from ACS imaging dScoville et all 
I2OO7I ; IScarlata et all l2007ah . we investigate the optical prop- 
erties of the galaxy population out to z ~ L Applying three 
straightforward classifications (spectral, photometric, and mor- 
phological), we are able to distinguish two main galaxy types, 
which are intimately connected with the bimodality and forma- 
tion history of the galaxy population: the three parametric clas- 
sifications constitute the axes of a "classification cube". It is im- 
portant to remark that we develop our cube using only observa- 
tional parameters, such as emission-line equivalent widths, con- 
tinuum indices, photometric colours (and not synthetic colours), 
and a non-parametric morphological classification. The choice is 
motivated by the desire to avoid dependence on models, galaxy 
templates, and parametrization of galaxy properties. Moreover, 
using simple observed properties makes the classification cube 
easy to replicate in other galaxy surveys. 

The principal aim of this paper is the definition of the galaxy 
classification cube, using a fair but relatively small sample of 
~ 1000 galaxies. The full exploitation of the classification cube, 
applied to larger galaxy samples (i.e. the lOK zCOSMOS sam- 
ple, Lilly et al. 2008, in preparation), will be presented in a forth- 
coming paper (Coppa et al., in preparation). 



2. Description of zCOSI\/IOS observations and data 
reduction 

For an exhaustive discussion of the zCOSMOS survey design, 
in particular of the assembly of the input catalogues, the tar- 
get selection, and the spectroscopic mas k design, w e address the 
reader to the survey presentation paper dLillv et all 12007). 

The brighter part of zCOSMOS project consists of spectro- 
scopic data of targets identified by means of a pure magnitude- 
limite d selection Iab < 22.5, as used in the CFRS (Lilly et all 
|1995 >) and VVDS-wide surveys ( Le Fevre et"an.l2005h . This se- 
lection culls galaxies mainly in the redshift range 0.1 < z < 1.2. 
A medium resolution grism (R ~ 600) has been used with a slit 
width of 1 arcsec, to achieve a velocity accuracy of ~ 100 kmi"^ 
and with one hour integrations to enable redshifts to be measured 
with a high success rate. The spectral range of observations is 
typically 5550 - 9650 A, which enables important spectroscopic 
diagnostic lines, centered on a rest-frame wavelength of approx- 
imately 4000A, to be measured out to a redshift of about 1 . The 
one hour total integration is acquired in five 720 s exposures, 
between which the telescope was placed at different offset po- 
sitions to complete an optimal sky subtraction, in a pattern of 
positions separated by one arcsec along the slit. The observa- 
tions were executed during periods when the seeing was better 
than 1.2 arcsec. The first zCOSMOS-bright spectroscopic obser- 
vations were completed in VLT Service Mode between April and 
June 2005. Eight sets of VIMOS mask observations were com- 
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pleted, which provided spectra of 1303 objects, corresponding to 
about 6% of the planned bright survey program. The data reduc- 
tion was carried out using the VIPGI software, developed specif- 
ically f or data acquired by VIMOS instrument ( Scodeggio et al.L 
l2005h . iLillv et all (12007) described the redshift determination 
procedures and the assignment of confidence classes to each 
measured redshift. The confidence in the redshift determination 
was indicated by a flag, whose values were: 4 (completely se- 
cure redshift); 3 (very secure redshift, but with a very marginal 
possibility of error); 2 (a likely redshift, but with a significant 
possibility of error); 1 (possible redshift); (no redshift deter- 
mination); and 9 (redshift based on a single strong emission 
line). The high quality of the photometric measurements for the 
COSMO S field and the correspo nding photometric redshifts at 
z < 1 .2 (iFeldmann et al.L l2006l) . implied that the photometric 
redshifts could be used to assess the spectroscopic redshift relia- 
bility of our zCOSMOS-bright spectra. To indicate the reliability 
of a redshift, we added a decimal place to its confidence flag: if 
the diff'erence between the photometric and spectroscopic red- 
shift, \zp - ZsV was less than 0.1 x (1 H- Zs), we added the decimal 
place=.5; otherwise, we added a decimal place=.l. In the fol- 
lowing, we include in the spectroscopic "high quality" data set 
all the galaxy spectra with flag 3 and 4, and those with flag 2.5, 
which includes galaxy spectra with flag=2 confidence for which 
the photometric and spectroscopic redshifts are in good agree- 
ment. We exclude from the high quality data set galaxy spectra 
of flag 9, where the redshift could be determined using only one 
emission line; although this flag corresponds to significant de- 
tection of an emission line, the absence of other spectroscopic 
diagnostic lines implies that for many the continuum signal-to- 
noise ratio (hereafter S/N) could be rather low. 



3. Spectral Analysis 

3. 1 . Spectral Measurements 

We performed spectral measurements for a total of 11 30 galaxies 
(and 20 QSOs) observed during the first year of observations and 
for which a redshift had been determined. This set of extragalac- 
tic objects forms the analyzed Ik zCOSMOS-bright sample. 

We measured the spectral quantities using semi-automatic 
procedures that were adaptations of the IRAF task splot. First, 
the continuum was fitted automatically to fixed wavelength in- 
tervals, although our procedures also enabled interactive adjust- 
ments of the continuum level to improve the line measurement 
in noisy spectra. Then, equivalent widths (EWs) and fluxes were 
measured by both applying a Gaussian-function fitting algorithm 
and direct integration of the continuum-subtracted line profiles. 
The fluxes measured by both techniques were in excellent agree- 
ment: for all emission lines, the two measurements agreed to 
within 25%, while the l-cr dispersion of the difference distribu- 
tion corresponded to only 11%. Figure [T] shows the distributions 
for principal spectroscopic line measurements for the lA:-sample. 
The equivalent widths are given in the rest-frame, and indicated 
by positive or negative values for emission or absorption lines, 
respectively. Errors in the measured equivalent widths were es- 
timated by taking into account both the measured rm.s. of the 
continuum close to the lines and the difference between the val- 
ues derived by the two techniques. The median uncertainty in 
an individual EW is 0.7 A, but it is worth noting that we did 
not account for the stellar absorption in the measurements of 
the Balmer emission lines, since in most cases the S/N of our 
spectra did not enable a simultaneous fit of both the absorp- 
tion and emission components. We measured EWs and fluxes 
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Fig. 1. Rest-frame measurement distributions of the analyzed 
spectral features, from top-left to bottom-right: D4000, and 
EW of the Unes [OII]/l3727, [OIII]^5007, H5, H/3, Ha, and 
[NII]/}6584. In each panel, we indicate the number of galaxies 
for which the feature is within the spectral observed range and 
the number of measured lines. 

for [OII]/l3727, H6, li/3, [OIII]^5007, Ha, and [NI I]^6584 
asweU as the strength of the 4000 A break (D4000; 'Bruzuall 
Il983h . The continuum fluxes were measured in fixed wavelength 
ranges, and our procedure computed the average using a sigma- 
clipping technique, to ensure that spikes, due to bad sky subtrac- 
tion and/or to cosmic-ray residuals, did not affect significantly 
the measured values. Errors in the continuum fluxes were then 
computed from the standard deviation. We focused our analy- 
sis on EW measurements and flux ratios because the absolute 
spectrophotometry was not completely reliable due to the cor- 
rections for slit losses that had been applied to the spectra. This 
set of "semi-automatic" measurements was used, in addition, to 
test and calibrate the ful ly automatic measurement code PlateFit 
jLamareille et al.Ll2006l) . which is being applied to larger zCOS- 
MOS data sets. 

In order to estimate our emission line detection limits, we 
measured the signal-to-noise ratio in the continuum adjacent to 
the lines. Figure |2] shows the relation between the S/N in the 
continuum and the observed equivalent widths of the measured 
emission lines in our sample: from this figure, it is clear that a 
correlation exists between the minimum detected EW and the 
continuum S/N. The lower envelope of the measured EWs (dot- 
ted line in Fig.|2]i is well described by the curve 



EW(detection limit) 



SL*A 
iS/NU„ 



(1) 



where A is the resolution element (in A) of our spectra and 
SL is the significance level of the detect able line, expressed 
in terms of sigma of the continuum noise (IColless et al.L Il990l : 
Manning, 2002). On the basis of Fig.|2] we adopt aSL(a) = 3.5. 
Measuring the S/N in the continuum of each galaxy spectra with- 
out a detectable emission line, we estimate the EW upper limits 
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Fig. 2. Observed EWs of all the measured emission lines (see 
Fig.[T]| versus the signal-to-noise ratio of the continuum close to 
the line; the lower envelope of the measured equivalent widths is 
fitted well by the dotted curve (see Eq.[T]in the text). The larger 
filled square symbols represent the [Oil] EW values measured 
in composite spectra of galaxies that individually did not demon- 
strate a detectable emission line, grouped in four equally popu- 
lated bins according to their continuum S/N; the line EWs have 
been measured in the composite spectra shifted to mean redshift 
of the corresponding galaxy group. The horizontal error bars in- 
dicate the range of S/N included in each bin, whereas the vertical 
error bars show the variation in the observed EW if the compos- 
ite spectra were moved to the extremes of the redshift range. 



(using Eq.[T]l. The robustness of this technique in estimating the 
line upper limits is demonstrated by the detection of the emission 
lines (in particular [OII]/l3727) in composite spectra of galax- 
ies without a detectable emission line. We grouped 141 single 
spectra in four equally populated bins according to their contin- 
uum S/N, and all the measured [Oil] EWs lie below the curve 
calculated using the above equation (large squares in Fig.|2|. 



3.2. Galaxy Spectral Classification 



Following iMignoli et al] (l2005h . we appUed to the zCOSMOS- 
bright galaxy data set the straightfo rward spectral classification 
scheme adopted for the 7^20 survey jCimatti et al.L 12002). which 
uses the equivalent width of the [OII]/i3727 line and the 4000 A 
break index to define two main spectral classes: red quiescent 
galaxies, which show a large D4000 value and a faint, if not 
undetected, [Oil] emission line; and blue star-forming galax- 
ies, with small D4000 and intense emission lines. This spectral 
classification scheme naturally arises from the intrinsic bimodal 
nature of galaxy physical p roperties that has been observed 
in the colour distribution dBaloghetall l2004bt iBelletal ' 



2004') and mean spectroscopic indices (^Kauff mann et al.L 12003 



f ranzetti et al., 2007). In Fig.[3j the distribution of the rest-frame 
[Oil] equivalent widths is shown for the zCOSMOS galaxy lA:- 
sample. A clear bimodality is evident in the EWo([OII]) dis- 



Fig.3. Distribution of the rest-frame [Oil] EWs for the zCOS- 
MOS galaxy lA:-sample. The gray histogram represents the pa- 
rameter distribution for galaxies included in the high-quality 
data set, with leftward arrows indicating the upper limits cal- 
culated using Eq. [T] the thin empty histogram is for the whole 
sample. 



tribution, especially if we consider that the histogram with left- 
ward arrows represents the upper limits to line measurements. 
On the basis of this behaviour, we separate the galaxy popula- 
tion into two classes, strong and weak line emitters, divided by 
an EW()([OII])=5A. With this threshold, almost all the objects 
with undetected emission lines are included in the weak emis- 
sion line galaxy group. Most galaxies with high EW that are not 
included in the high-quality sample are flag=9 objects (see end 
of Sect. 2). 

Looking at Fig.|4] a couple of comments can be made. First, 
the "low-density" region in the leftward corner of the classifica- 
tion plane is populated mainly by poor-quality flagged objects 
(indicated by crosses in the figure), confirming both the doubt- 
ful nature of their redshift identification and the efficiency of 
the spectral classification scheme. The galaxies in this region 
of the plane are not spectroscopically classified. Furthermore, 
when applied to the high-ranked zCOSMOS galaxies, the clas- 
sification scheme demonstrates its efficiency: only four objects 
fall marginally in the leftward region of the plot, whereas the 
intermediate class (red emission line galaxies) constitutes only 
4% of the sample. Interestingly, the few (broad line) active nu- 
clei observed in the useful redshift range are quite segregated 
in the classification plane (see the encircled symbols in Fig.|4]i, 
having very blue continuum but relatively faint [Oil] emission 
lines (<10A, Forster et al.Ll200ll) . 

Figure |4] shows the resulting spectral classification plane. 
Due to the high quality of zCOSMOS spectroscopic data, higher 
than in the K2Q survey, we are able to appreciate how most of the 
galaxies with intense emission lines occupy a narrow region in 
the D4000-[OII] plane. We performed an iterative linear least 
squares fit with cr-clipping, to define and delimitate the locus 
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Fig. 4. The classification plane for the Ik zCOSMOS galaxies; using the rest-frame [OII]/l3727 EW and D4000, we define three 
main classes: red quiescent galaxies (lower right region), blue emission-line galaxies (upper left region) and the intermediate 
galaxies, which show emission lines but a large value of D4000. In the figure, we plot all classifiable galaxies (with z ~ 0.45 -Hi. 25), 
with different symbols following the redshift quality ranks (crosses for the low-quality spectra, filled squares for the high-quality 
ones). Broad line objects (quasars) are encircled in figure. The dashed lines represent the boundaries of the adopted classification 
scheme following Eq. |2] The percentage of high quahty galaxies falling in each classification regions are labelled at plot corners. 
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Fig. 5. Composite spectra of the three main galaxy classes as defined in Fig. H] From bottom to top: 126 red quiescent, 20 interme- 
diate and 361 blue star-forming galaxies contribute to the averages. 
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of these objects (see diagonal dashed lines in Fig. |4]i. The final 
criterion adopted to define the blue emission line galaxies was 

EWo([011]) > 5k, 
1.50 < D4000 + 0.33 x log(£Wo([OII])) < 2.22 (2) 

A relationship between the D4000 and the [Oil] equivalent 
width was expected, since both features are related to the galaxy 
star-formation history, either the most recent one (£'Wo([OII])) 
or its integral over time (D4000). It is particularly encouraging 
that, by projecting our linear fit in the D4000-specific SFR plane, 
with the simple assumption of constant mass-to-luminosity ratio, 
we reproduce almost exactly the mode of the SPSS da ta distri- 
bution presented in Fig. 1 1 of iBrinchmann et aTl (1200 4V 

It is worth keeping in mind that the observed spectral range 
enables us to apply this spectral classification scheme only to 
galaxies of sufliciently high redshift to include the [OII]-Z)4000 
region, i.e. the redshift range z ~ 0.45 - 1 .25. Within this redshift 
range, 63 1 galaxies were classified, of which 511 belong to the 
high-quality sample on the basis of the assigned redshift flag. 
More than 80% of the spectroscopically classified zCOSMOS 
galaxies belong to the two previously defined classes: red qui- 
escent galaxies and blue star-forming galaxies. A third group 
of objects, which have emission lines but a red stellar contin- 
uum (as implied by their high value of D4000), constitutes what 
we refer to as the intermediate class of galaxies. Interestingly, 
20% of these objects show some sign of nuclear activity: three 
of the twenty galaxies in this class are classified as Seyfert 2 
galaxies (Bongiorno et al., in prepara tion) on the basis of the di- 
agnostic diagrams of Baldwin et al.l (11981). and on e additional 
galaxy has an X-ray counterpart dBrusa et al.Ll2007h . In contrast, 
the fraction of the blue star-forming galaxy population that are 
candidate AGN is only 2%, i.e. 7 out of 361 galaxies. 

Composite spectra for each galaxy class were generated by 
averaging all spectra included in that class. To create the com- 
posites, each spectrum was shifted to the rest-frame according 
to its redshift (with a 1 .5 A rest-frame birQ) and normalized in 
the 4000-4500 A wavelength range, which was always present 
in the observed spectroscopic window. An identical weight was 
assigned to each individual spectrum, to avoid biasing the fi- 
nal composite towards the brightest galaxies. Finally, all spec- 
tra belonging to a particular spectroscopic class were stacked. 
In Fig. |5] the average spectra of the three main spectroscopic 
classes are plotted. The continuum in the composite of the emis- 
sion line galaxies with red D4000 is largely indistinguishable 
from that of the purely passive galaxies composite. Therefore, 
this class of intermediate objects should consist mainly of galax- 
ies undergoing a modest star-formation episode, rather than of 
heavily reddened star-forming galaxies. The presence of the 
[OIII] emission doublet, with no evidence of Hp in emission, in 
the composite of intermediate objects could indicate some nu- 
clear activity (i.e. type II AGNs) con tamination in th is galaxy 
class, as already noted above (see also lYan et al.Ll20()6h . 

4. Photometric colours of galaxy spectral classes 

The bimodality in the galaxy spectroscopic index measurements 
identified above may be naturally linked to the colour bimodality 
observed widely in the galaxy population. We therefore compare 

' at the median redshift of the analyzed galaxy sample (z = 0.724), a 
rest-frame bin of 1.5 A matches the pixel size used in the spectroscopic 
observations. 



the observed colours and spectral classification of galaxies that 
belong to the zCOSMOS-bright U-sample. 

The photometric observations of the COSMOS field that we 
use in our study include optical (Bj, g^ , Vj, r^, and data 
acquired using SuprimeCam at the Subaru Telescope, M*-band 
data with the Canada-France-Hawaii Telescope and /814 images 
with ACS/HST) and near-Infrared (/Tj-band data acquired us- 
ing the Cerro Tololo International Observatory and Kitt Peak 
National Observatory telescopes) wavelengths. Details of the 
ground-based observati ons and data reduction are presented in 
Capak et al. (2007) and Taniguchi et an(l2007h . The reduced im- 
ages in all bands were PSF-matched by Gaussian convolution 
with a final FWHM corresponding to the seeing in the K^- 
band (1.5"); the multi- band photometry catalog was then gener- 
ated using SExtractor (i Bertin & Arnoutsl Il996l) . First, the total 
(SExtractor magauto) and fixed-aperture (3"-diameter) magni- 
tudes were measured on the detection image (i-band), allowing 
the estimate of the aperture correction for each galaxy. This cor- 
rection was applied subsequently to aperture magnitudes mea- 
sured in the other optical/near-IR bands. Details of the photom- 
etry, star/g alax y separ ation, and catalog construction were pro- 
vided bv lCapak etal.i (i2007b . 

We analyzed diff'erent colours from the photometric catalog 
and attempted to identify a separation in photometric space be- 
tween the above defined spectral galaxy classes. The B-z colour, 
in addition to both B - K and u - i, appears to provide the most 
effective way of separating red quiescent and blue star-forming 
galaxies. In Fig. |6] we illustrate, in the B-z distribution his- 
tograms (upper panel; cyan histogram for blue star-forming, red 
histogram for quiescent galaxies, both with and without emis- 
sion lines), the good separation between the colour distribution 
of the different spectral classes, although with some overlap for 
2 < B-z < 3. We verified whether the slit aperture effects would 
bias the comparison between the observed colours and spectral 
classification: the slit corrections, estimated on the basis of the 
difference between the integrated magnitudes computed from the 
flux-calibrated spectra and the I-band photometric data, do not 
show any trend both with the magnitude and B-z colour in the 
analyzed sample. 

The correlation between colour and spectral classification 
becomes even tighter when we consider galaxy redshifts. The 
B-z versus redshift plot (lower panel; cyan and red squares 
for star-forming and quiescent galaxies, respectively, magenta 
crosses for intermediate ones) exhibits a clear segregation be- 
tween the galaxy classes. To divide eff'ectively the two galaxy 
populations, we adopted the B-z versus redshift evolu- 
tionary track of a template spectrum from the Coleman, Wu 
& Weedman Extended set (IColeman. Wu & WeedriianL 119801: 
Sawicki et al... 1997) with properties corresponding most closely 
to those dividing th e two populations. The chosen template had 
the o ptical col ours ("Fukugi ta et al.L 1 19951) and spectral appear- 
ance dKennicu tt, 1992) of an Sab galaxy. Also in local galaxy 
samples (i.e. T he Revised Shapley-Ames Catalog of Galaxies; 
Ivan den Berglu , .200 7), the transition between blue and red galax- 
ies occurs at the Hubble type Sab/Sb. 

The lower panel of Fig. |6]shows that almost all of the small 
number of discordant galaxies (blue quiescent and red star- 
forming) fall close to the separation track in the colour classi- 
fication plane, confirming the agreement between the spectral 
and photometric classification. Moreover, galaxies of the inter- 
mediate spectral class (red D4000 with emission lines, magenta 
crosses in Fig.|6|l have a colour distribution similar to that of qui- 
escent galaxies, in agreement with the observed similarity of the 
stellar continuum in the composite spectra noted above (Fig.|5]l. 
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Fig. 6. Upper Panel: observed B-z colour distribution of the galaxies belonging to zCOSMOS-bright 1^-sample with 0.45 < z < 1 .2; 
cyan histogram for blue star-forming, red histogram for quiescent galaxies, both with and without emission lines. Lower Panel: B-z 
versus redshift plot for the whole high quality sample; cyan and red squares for star-forming and quiescent galaxies respectively, 
magenta crosses for intermediate ones). We divide the two galaxy populations using as discriminator the B-z track of the Sab 
template from the Coleman Wu & Weedman Extended set (solid line in figure: , Sawicki et al..(.1997.) ). 



By merging these two spectral classes and using the track 
of the Sab-template to separate the photometrically "blue" and 
"red" galaxies, we compiled a 2 x 2 contingency table for the 
high quality sample (Table [U, useful in performing a statistical 
analysis of the relationship between the spectral and photomet- 
ric classification. To quantify the degree of correlation between 
the two classification schemes, we computed Pearson's linear 
correlation coefficient for the data presented in Table [T] the two 
classifications are well correlated, as indicated by a Pearson co- 
efficient of 0.77. 

Table 1. Numbers of galaxies in the spectroscopic and photo- 
metric classes 





Spectral classification 


B-z 


Quiescent Star-forming TOTAL 


Red 


119 20 139 


Blue 


27 335 362 


TOTAL 


146 355 501 



In the following, we summarize how the different galaxy 
populations inhabit the colour versus spectral classification 
plane: 



- 82% of objects classified as "quiescent" from spectra are red; 

- 94% of objects classified as "star-forming" from spectra are 
blue; 

- 86% of objects with red colours are "quiescent" from the 
spectra; 

- 93% of objects with blue colours are "star-forming" from the 
spectra. 

- The fraction of outHers is -9% (47/501). 

The correspondence between the two classifications is not opti- 
mal for the 18% quiescent galaxies that show bluish colours. We 
investigated the nature of these objects by studying their com- 
posite spectrum. The continuum in the composite of blue qui- 
escent galaxies is largely indistinguishable from that of the red, 
quiescent galaxies, but the spectrum shows a moderate [Oil] 
emission line {EWo = 4.4A) and a strong H6 in absorption 
(EWo = -3.8A). The presence of these two spectral features 
represents star formation activity, which could also explain the 
bluish colours of this class of galaxy. The rest-frame equivalent 
width of [Oil], measured in the composite spectrum, is below 
the cut adopted for the spectral classification (£'VVo([OII]) > 
5A), and only one of the 27 blue quiescent galaxies is spectro- 
scopically classified as intermediate. Therefore, it is more rea- 
sonable to consider these objects as being galaxies in-between 
the two main spectral classes, rather than galaxies of extremely 
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discordant spectral and colour properties. Indeed, about 30% 
of these galaxies can be spectroscopically classified as post- 
starburst galaxies and constitute the 64% of the Ik sample of 
post-starburst galaxies (7/11; Vergani et al., in preparation). 

The good agreement between the classification from spectral 
data (quiescent/star-forming galaxies) and that based on colours 
(red/blue galaxies) is achieved within the relatively bright (Iab < 
22.5) z COSMOS samp le; using the deeper {Iab < 24.0) VVDS 
survey. IFranzetti et alj (12007 ) demonstrated that the incidence of 
red star-forming or blue quiescent galaxies appear to increase 
when studying fainter galaxy samples. 



5. ACS Morphology of galaxy spectral classes 

The HST/ACS imaging is an integral component of the 
COSMOS survey, providing very high sensitivity and high res- 
olution imaging over a large field of view. These space-based 
images yield resolved morphologies for several hundreds of 
thousands galaxies, and the zCOSMOS-bright survey magnitude 
limit (Iab < 22.5) is matched with the I-band ACS depth to allow 
classical bulge-disk decomposition. The COSMOS survey is a 
natural "battlefield" for competing morphological classification 
schemes. Among all the possible and equally good techniques 
implemented by the various teams within the COSMOS col- 
laboration, we selected as our morphological class ification tool 
the Zu rich Estimator of Structural Types (ZEST; IScarlata et al.L 
l2007a[) . We adopted this approach mainly because ZES T classi- 
fies galaxy types by applying a variety of non-parametric tech- 
niques that measure physical parameters directly on the ACS im- 
ages in a way similar to how we classify galaxies in terms of both 
spectroscopic and optical colour data. 

In brief, ZES T quantitatively describes the galaxy structure 
using three variables (PCI; PC2; PC3), obtained by performing 
a principal component analysis (PCA) in the five-dimensional 
parameter space defined by asymmetry (A), concentration (C), 
Gini coefficient (G), the 2nd-order moment of the brightest 20% 
of galaxy pixels (M20), and the ellipticity of the light distribution 
(e). The morphological classification is performed in the PCl- 
PC2-PC3 space by linking a (dominant) morphological class to 
different regions in the PC-space. Finally, ZES T assigns each 
analyzed object to a straightforward morphological Type {- 1 
for early-type galaxies, = 2 for disk galaxies and = 3 for ir- 
regular galaxies). Furthermore, a "bulgeness" parameter (from 
for bulge-dominated galaxies to 3 for bulge-les s disks) is also 
assign ed to each galaxy classified as Type 2 (IScarlata et al.L 
l2007ah . following the median value of the Sersic indices distri- 
bution of galaxies clustered in the PC-space. The Sersic indices 
were computed by Sargent et al. (2007). 

We compare the spectral classification with the morphologi- 
cal characterization assigned by ZEST in Fig.|7] where the his- 
tograms represent the frequency distributions of the morpholog- 
ical types for the three main spectroscopic classes. This figure 
illustrates a fairly good correlation between the spectral and the 
morphological classification, represented by a gradual increase 
in the fraction of star-forming galaxies towards morphological 
late-types, and, similarly, an increase in spectroscopic quiescent 
galaxies towards the morphological early- types. 

A more instructive way to look at the correspondence be- 
tween the morphological and spectral classifications is pre- 
sented in Fig. [8] which shows the distribution of the mor- 
phological types in our spectral classification plane. Following 
the morphological characterization, we divided the galaxies 
into: (/ ) "spheroids", a class including both the pure ellipticals 



Biilse Domin. Pure Disk 



cr 




Ellipiical 



Disky 
ZEST_type 



rr'cgular 



Fig. 7. Frequency distributions of ZEST morphological types 
for the spectroscopic blue emission line galaxies (top his- 
togram), the spectroscopic intermediate galaxies (i.e. galaxies 
which show red continuum with emission lines; middle his- 
togram), and the quiescent galaxies (bottom histogram). The 
four central bins represent the most-populated ZES T class of 
disk galaxies, further splitted on the basis of the "bulgeness" pa- 
rameter The numbers on top of the bars indicate the number of 
galaxies in each bin. 



{ZEST -type = 1) and the bulge-dominated galaxies (ZEST- 
type = 2 with bulgeness parameter = Oj3; (H) "morphologically- 
mixed" objects, disky (ZES T-type - 2) galaxies with interme- 
diate bulgeness parameter {- 1); (///) "disk-dominated" galax- 
ies with ZEST-type = 2 and bulgeness parameter > 2; and 
[iv ) "irregulars", with Z/s^ T-type = 3. The four panels of Fig.|8] 
demonstrate that all but one of these morphologically defined 
galaxy types appear to occupy the regions expected based on 
their spectral classification. More than 90% of the morphological 
ellipticals are quiescent from the spectra, as well as a large frac- 
tion (82%) of the spheroids. Also, more than 90% of the disk- 
dominated and irregular galaxies have blue emission line spec- 
tra. Only the morphologically mixed objects (i.e. disk galaxies 
with intermediate B/T parameters) are spread all over the spec- 
tral classification plane. Table|2]shows the distribution of ZES T 
types in the lA:-galaxy sample as a function of their spectroscopic 
class. 

The application of a 2 x 2 contingency table, as a statisti- 
cal test of correspondence between the spectral and morpholog- 
ical classifications, is not a straightforward exercise in this case, 
since it is difficult to assign objects with ZEST-type = 2.1 to an 
ear ly- or late- galaxy type. 

iBailin & HarrisI (l2008l) discovered a similar population of 
morphologically intermediate SDSS galaxies, finding that in- 
clination effects play an important role in their classification. 
Nevertheless, both a visual inspection of the ACS snapshots and 
a statistical analysis of the asymmetry and ellipticity parame- 



- the inclusion of the bulge-dominated disk galaxies into the same 
class as pure ellipticals was already proposed bv lScarlata et alJ (l2007bl) 
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Fig. 8. Distributions of different morphological groups in our 
spectral classification plane. Panel a): Spheroidal galaxies, in- 
cluding both pure ellipticals {ZES T-type = 1 ; filled circles) and 
bulge-dominated galaxies (ZES T-type - 2.0; empty squares). 
Panel b): disky galaxies with intermediate bulgeness parameter 
(ZE^T-type =2.1; empty squares). Panel c): "disk-dominated" 
galaxies {ZES T-type = 2.2 and = 2.3; empty squares). Panel d): 
Irregular galaxies (Z/s^r-type = 3; filled triangles). In all the 
four panels, the dashed lines indicate the spectral class bound- 
aries. 

Table 2. Numbers of galaxies in the spectroscopic and ZEST 
classes 





Spectral classification 


XES T-type 


Quiescent Star-forming TOT\L 


1.0-2.0 


76 24 100 


2.1 


55 70 125 


2.2-3.0 


1 1 256 267 


TOTAL 


142 350 492 



ters do not suggest the presence of an incUnation bias in the 
zCOSMOS galaxies with Z/i^T-type =2.1. This intermediate 
morphological class is plausibly composed of galaxies spanning 
a wide range of bulge-to-disk ratios and ZES T alone is unable 
to identify a dominant component in their visual morphology. 
We decided to fully utilize the classification cube potentialities, 
including the photometric colours, in the analysis of this inter- 
mediate morphological class. Figure |9] shows the colour dis- 
tribution of Zii^r-type =2.1 galaxies, divided following the 
spectral classification. The upper panel shows the distribution 
of the differences between the observed {B - z) colour and the 
colour, at the same redshift of each galaxy, of the Sab template 
adopted to separate the blue and red populations in the pre- 
vious section, for galaxies spectroscopically classified as star- 
forming; the lower panel shows the same colour difference for 
the Type = 2.1 galaxies spectroscopically classified as quiescent. 
The figure clearly shows that galaxies with ZEST -type =2.1 
have B - z colours largely in agreement with their spectral clas- 



Fig. 9. The distribution of the differences between the observed 
(B - Z) colour of the Type =2.1 galaxies and the colour, at 
the redshift of each galaxy, of the Sab template adopted to pho- 
tometrically separate the blue and red populations in Fig. |6] In 
the upper panel, the galaxies spectroscopically classified as star- 
forming; in the lower panel, the galaxies spectroscopically clas- 
sified as quiescent. 



sification: 93% (65/70) of the ZEST Type =2.1 classified as 
"star-forming" from the spectra have blue colours, whereas 71% 
(40/55) of the ZEST Type = 2.1 classified as "quiescent" from 
the spectra have red colours. Interestingly, the latter galaxy class 
exhibits a bimodal colour distribution, which is separated well 
by the colour track of the Sab template. We investigated the aver- 
age properties of the two classes of red and blue quiescent galax- 
ies with intermediate ZEST morphology by means of the anal- 
ysis of their composite spectra: the stellar continuum appears 
similar in the two composites, resembling that of the normal qui- 
escent galaxies (see Fig.|5]l, even if the average spectrum of the 
blue quiescent class is of low signal-to-noise ratio because of the 
small number of coadded galaxy spectra. The main difference is 
the presence of emission lines in the blue quiescent composite 
(with EW([OII])~5A), while emission lines are undetected in 
the red quiescent composite. However, the careful analysis of 
single galaxy sub-classes is beyond the scope of this paper, and 
are the subject of a forthcoming paper in which the properties of 
these galaxy classes will be revisited using a larger zCOSMOS 
data set. 



Table 3. Numbers of galaxies in the spectroscopic and morpho- 
logical classes 



Morph+(B - z) 

Spheroids 
Disk/Irr 
TOTAL 


Spectral classification 
Quiescent Star-forming TOTAL 


115 29 144 
27 321 348 
142 350 492 
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We therefore separated the galaxies of this morphological 
class by colour, by assigning the red-coloured Type 2.1 galaxies 
to a Spheroid class and the blue-coloured Type 2.1 galaxies to 
a Disk/Irregulars category. Following this recipe, we derived a 
2x2 contingency table (Table |3]l, which provided a clear indica- 
tion of the correlation between the two classifications (Pearson 
correlation coefficient^ 0.73). 



6. The galaxy classification cube 

It is evident from the previous sections that strong correla- 
tions are present between the spectral classification and both 
the colour and morphological properties of the zCOSMOS- 
bright galaxies. In this section, we merge these classifications 
into a three-dimensional characterization, the galaxy classi- 
fication cube. We include in our analysis all galaxies with 
0.45 < z < 1 .25 Q and complete photometric and morphological 
information. The full sample includes 577 galaxies, of which 
487 (84%) belong to the high-quality spectroscopic sample on 
the basis of the assigned redshift flag. 

A three-digit code was then assigned to each galaxy of the 
sample, such that each digit represented one of the three classifi- 
cations (spectral, photometric, and morphological, respectively). 
Each digit can assume a bitwise value. This corresponds to 1, 
if the galaxy was spectroscopically classified as quiescent, had 
red B - z colour, or its morphology was dominated by the bulge 
component. Alternatively, we assigned a digit value of 2 if the 
galaxy was spectroscopically star-forming, had blue B-z colour 
or showed a disky or irregular morphology. 

We recall two caveats related to both the procedure we 
adopted to separate galaxies into one of the two classes, 
and the presence of galaxies that deviate from the bimodal 
trends in properties (quiescent/star-forming, red/blue, bulge- 
dominated/disk-dominated). First, a galaxy spectrum with emis- 
sion lines but a red stellar continuum (implied by a rela- 
tively high value of D4000) was classified spectroscopically 
as quiescent because their global properties indicated that they 
were probably ellipticals experiencing a modest star-formation 
episode or including a faint/obscured AGN. Secondly, we 
used the B-z colour to determine a more robust separation 
of the galaxies included by ZEST in the intermediate class 
(Z£5r Type = 2.1). 



Table 4. The Classification Cube: galaxy counts and class per- 
centages 



class-code 


High-Q. Sample 


All Galaxies 


s-c-m 


N 


(%) 


N 


(%) 


2-2-2 


311 


(63.9) 


374 


(64.8) 


1-1-1 


105 


(21.6) 


113 


(19.6) 


1-2-2 


18 


(3.7) 


26 


(4.5) 


2-1-1 


13 


(2.7) 


17 


(2.9) 


2-1-2 


7 


(1.4) 


8 


(1.4) 


1-2-1 


9 


(1.8) 


11 


(1.9) 


2-2-1 


15 


(3.1) 


17 


(2.9) 


1-1-2 


9 


(1.8) 


11 


(1.9) 


TOTAL 


487 


(100) 


577 


(100) 



^ since the features adopted for spectral classification (£)4000 and 
EW([OII])) can be seen in the zCOSMOS-bright spectra only within 
this redshift range 



In Table [H numbers and percentages of objects belonging 
to each cell of the classification cube are presented, using the 
three-digit code to identify the galaxy groups. The numbers for 
both the high-quality and the full sample are listed, and although 
all the single classification definitions were determined by using 
the high-quality sample, the percentages of galaxies included in 
each group vary only slightly between the two different sam- 
ples. The main result emerging from the 3D analysis of the op- 
tical properties of the zCOSMOS-bright galaxies is that a large 
fraction of them shows a fully concordant classification: 85.5% 
(84.4%) of the analyzed galaxies in the high-quality (full) sam- 
ple can be assigned consistently to the red or blue branch of 
the bimodal distribution using both a spectral, photometric, and 
morphological criterion. Due to the redshift constraints imposed 
by the spectral classification, we cannot explore a very large red- 
shift domain. However, we divided the useful redshift range into 
two bins, so that each bin contained a similar number of ob- 
jects: in the low-z range (0.45 < z < 0.73) 83.4% of the galax- 
ies shows a fully concordant classification, while in the high- 
redshift range (0.73 < z < 1.25) this percentage slightly grows 
to 85.5%. We conclude that the fraction of galaxies for which 
the three classification schemes are in agreement remains con- 
stant (within our statistical errors) for the entire redshift range 
covered by our analysis. 

7. Conclusion 

We have presented a classification scheme for galaxies observed 
by the Ik zCOSMOS-bright sample, developed using spectral, 
colour, and morphological information. Since the spectral clas- 
sification that we adopted is based on measurements of the 
[Oil] and D4000 features, and given the wavelength range of 
the zCOSMOS spectra, we applied this classification scheme to 
galaxies in the redshift range 0.45 < z < 1.25. We demonstrated 
that the adopted spectral classification is an effective way of sep- 
arating galaxies into two main classes, due to the high quality of 
the zCOSMOS data. 

We have found that an excellent correlation exists between 
the classification from spectral data (quiescent/star-forming 
galaxies) and that based on colours (red/blue galaxies), and that a 
similarly good correlation exists between the classification based 
on spectral data and that based on morphological analysis (early- 
/late-type galaxies) only if we complement the morphological 
classification with additional colour information. The "outliers" 
in this simple 2-way analysis of different classification meth- 
ods (i.e. spectra versus colours and spectra versus morphology) 
are quiescent galaxies (high D4000 value) with emission lines 
(about 4% of the total sample of galaxies) and morphologically 
intermediate galaxies (about 20% of the total sample of galaxies) 
that consists of both blue/star-forming and red/quiescent galax- 
ies. These are the morphologically-classified galaxies that re- 
quire additional colour information to ensure good agreement 
between the classification of schemes using both spectroscopic 
and morphological data. 

Finally, analyzing the distribution of all galaxies in a 3D- 
classification cube, based on the three different classification 
methods, we have found that for about 85% of the galaxies there 
is complete agreement between the classifications (i.e. quies- 
cent, red, ellipticals or bulge-dominated galaxies (~ 20%) or 
star-forming, blue, disk-dominated galaxies (~ 65%)). If addi- 
tional colour information was not used for the ZES T-type = 2.1 
galaxies (i.e. by assigning all of them to the spheroid- or to the 
disk-galaxies), the percentage of fully concordant classifications 
would decrease to about 75%. 
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Fig. 10. Classification cube and rest-frame U - V colour of the 
Ik zCOSMOS-bright data set. The upper panel shows the rest- 
frame colour distributions of galaxies with fully concordant clas- 
sification, the blue histogram representing the star-forming, blue, 
disk-dominated galaxies, whilst the red histogram is for the qui- 
escent, red, ellipticals or bulge-dominated galaxies. The lower 
panel shows the rest-frame colour distribution for the outliers, 
i.e. galaxies which do not present the same behaviour in all the 
three classifications. 

The correspondence between the classification cube and 
classical representation of galaxy bimodality, which uses the 
rest-frame U - V colour^ is shown in Fig. [TO] In the upper 
panel, as expected, the rest-frame U - V colour distributions of 
the two galaxy groups with fully concordant classifications are 
separated well and reproduce the bimodality. Perhaps more in- 
teresting is the distribution of outliers in the lower panel, which 
is centered in-between the two branches of the galaxy bimodal- 
ity, in the so-called 'green valley', the relatively sparse region 
between the blue cloud of star-forming galaxies, and the red se- 
quence of quiescent ones. 

The results presented in this paper demonstrate the poten- 
tial of the data set provided by the COSMOS/zCOSMOS sur- 
vey for investigating the nature of the observed relationships be- 
tween galaxy properties and highlight fundamental aspects of 
the galaxy population. The main conclusion of this work is that 
the bimodality is a prevalent feature of the galaxy population, 
at least out to z ~ 1, in three critical parameters describing 
galaxy physical properties: morphological, colour, and spectro- 
scopic diagnostic line data. More detailed analysis of the prop- 
erties of both a wider sample of galaxies and outliers in this 
3D-classification cube, based on an analysis of the larger lOK 
zCOSMOS-bright sample, will be presented elsewhere. 
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